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ABSTRACT 

Fjortoft»s  graphical  method  of  prognosis  of  the  500-mb 
surface  assumed  relative  vorticity  change  as  being  due  to 
advection  only.   This  study  considers  the  total  24-hour  rela- 
tive vorticity  change  as  consisting  of  two  components,  that 
due  to  100$  advection  and  that  due  to  non-advective  changes. 
Objective  methods  are  developed  for  the  prognosis  of  the 
latter.   These  changes  are  then  added  to  the  advective  changes 
to  give  the  total  24-hour  relative  vorticity  change  from  which 
the  24-hour  height  change  can  be  recovered. 

The  500-mb  velocity  field  is  subjectively  categorized 

into  five  types.  For  each  type  the  non-advective  relative 

vorticity  change  is  correlated  with  a  parameter  determined 

from  the  spreading  of  the  contours.   Results  for  the  five 

categories,  when  applied  to  data  for  the  months  of  February 

and  March  1957,  are  as  follows: 

Correlation 
Case  Velocity  Field  Coefficient 

I      Marked  Cyclonic  Curvature  0.60 

II      Moderate  Cyclonic  Curvature  0.71 

i 

III      Negligible  Curvature  0.79 

IV      Moderate  Anticyclonic  Curvature  0.56 

V      Marked  Anticyclonic  Curvature  0.10 

Quantitative  objective  forecast  rules  are  developed  from  the 
regression  line's  relating  the  non-advective  relative  vorticity 
change  and  the  parameters  for  each  of  the  five  cases. 

Using  the  developed  technique  the  results  of  five  24-hour 
500-mb  test  forecasts  for  the  United  States  are  presented. 
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These  tests  show  an  Improvement  of  four  feet  in  forecast 
height  value  for  the  entire  forecast  area  as  compared  to 
the  FJortoft  method.  However,  in  considering  only  those 
points  at  which  the  computation  was  made  an  average  improve- 
ment of  34  feet  is  noted* 

Finally  the  persistence  of  non-advective  relative 
vorticity  changes  is  considered  for  its  value  as  a  height- 
change  forecast  tool.   The  results  of  seven  forecasts  pre- 
pared using  this  assumption  indicate  an  improvement  of  32 
feet  over  the  Fjortoft  method  for  the  entire  forecast  area* 

The  writer  wishes  to  express  appreciation  to  Professor 
W.  D.  Duthie  for  guidance  and  assistance  in  the  preparation 
of  this  paper* 


ii  (a) 


TABLE  OF  CONTENTS 


Page 
i 


CERTIFICATE  OF  APPROVAL 

ABSTRACT  1± 

TABLE  OF  CONTENTS  Hi 

LIST  OF  ILLUSTRATIONS  iv 

TABLE  OF  SYMBOLS  AND  ABBREVIATIONS  v 

CHAPTER  I      INTRODUCTION  1 

CHAPTER  II     THEORETICAL  DEVELOPMENT  4 

CHAPTER  III    PROCEDURE  8 

CHAPTER  IV    FORECASTING  VORTICITY  CHANGE  10 

CHAPTER  V      INDEPENDENT  TEST  RESULTS  26 

CHAPTER  VI    FURTHER  STUDIES  28 

BIBLIOGRAPHY  31 


iii 


LIST  OP  ILLUSTRATIONS 

Figure  Page 

lfl   Contour  pattern  for  Case  I  and  method  of  12 
measuring  x^# 

2.  Regression  line  of  R  on  xj..  14 

3.  Regression  line  of  R  on  X£«  16 

4*   Contour  pattern  for  Case  III  and  method  17 
of  measuring  x3. 

5*   Regression  line  of  R  on  x^*  18 

6«   Contour  pattern  for  Case  IV  and  method  of  20 
measuring  X&* 

7*   Regression  line  of  R  on  X4*  21 

8«  Regression  line  of  R  on  X5.  '      22 

Table 

I.   Summary  of  results.  23 


It 


TABLE  OP  SYMBOLS  AND  ABBREVIATIONS 

J        Vertical  component  of  the  relative  vorticity* 

7    Coriolis  parameter. 

t       Time •  / 

V   Horizontal  vector  velocity  of  the  wind. 

CJ       Individual  rate  of  change  of  pressure  with  respect 
to  time* 

R  Non-advective  relative  vorticity  change. 

Q  Gravitational  force* 

<//  Grid  distance* 

£  Contour  height* 

£  Space-mean  contour  height* 

Va  Horizontal  vector  velocity  of  the  geos trophic  wind* 

1/   Horizontal  vector  velocity  of  the  space-mean 
j       geos trophic  wind* 

x^  Case  I  divergence  parameter* 

*2  Case  II  divergence  parameter* 

X3  Case  III  divergence  parameter* 

X4  Case  IV  divergence  parameter* 

X5  Case  V  divergence  parameter* 

0         » $ 

jsr&      Earth»s  rotational  velocity. 
pf)      Map  scale  constant* 
Latitude* 


C*      :     j-k    «    fW) 


CHAPTER  I 
INTRODUCTION 
Fjortof t  (^  lj  developed  a  graphical  method  of  fore- 
casting the  500-mb  constant  pressure  surface  by  the  follow- 
ing procedure:   1)  determine  the  distribution  of  relative 
vorticity  at  the  500-mb  level;  2)  advect  this  vorticity  with 
the  mean  500-mb  flow;  3)  obtain  the  24-hour  vorticity  change 
by  subtraction  of  these  two  fields;  4)  recover  the  24-hour 
height-change  field  from  the  24-hour  vorticity-ohange  field; 
and  5)  arrive  at  the  prognostic  chart  by  adding  the  height- 
change  field  to  the  current  500-mb  chart.   The  assumptions 
involved  in  this  procedure  include  barotropy,  the  geostrophic 
approximation,  the  neglect  of  the  change  of  coriolis  para- 
meter and  map  scale  with  change  in  latitude,  and  the  conser- 
vation of  the  mean  500-mb  flow  during  the  forecast  period. 

Many  less  restrictive  models  have  been  developed. 
Fjortof t,  in  the  paper  cited  above,  compensated  for  the  change 
of  coriolis  force  and  map  scale  with  latitude  by  adding  a 
field,  which  is  a  function  of  latitude  only,  to  the  relative 
vorticity  field  and  advecting  this  quantity.   In  another  paper, 
Fjortof t  I  2  I  has  developed  a  baroclinic  model  that  assumes  a 
functional  relationship  in  the  wind  profile.   Hesse  j_3 J  has 
included  baroclinity  with  vertical  motion  due  to  terrain. 
Without  question  these  more  sophisticated  graphical  tech- 
niques will  prove  superior  to  arbitrarily  selected  statisti- 
cal parameters.  However,  for  the  meteoroligist  in  the  field, 


the  time  involved  in  the  preparation  of  the  routine  fore- 
oast  must  be  held  within  strict  limits.   It  is  felt  that 
techniques  more  complicated  than  the  simple  Pj or toft  method 
would  not  be  within  the  capabilities  of  the  average  Naval 
Weather  Station. 

The  Air  Weather  Service  [4]  has  tested  the  simple 
Fjortof t  method  as  a  means  of  objectively  forecasting  the 
500-mb  contour  field  for  a  24-hour  period.   In  the  test 
consisting  of  541  samples  the  Pjortoft  method  gave  an 
average  absolute  error  of  208  feet  in  the  forecast  500-mb 
height,  whereas  subjective  techniques  employed  at  the  TJSAF 
Weather  Central  gave  an  average  error  of  207  feet.   It  was 
concluded  that: 

objective  use  of  the  graphical  method  presented 
by  R.  Pjortoft  produces  24-hour  500-mb  prognostic 
charts  of  about  the  same  quality  as  those  produced 
by  skilled  forecasters  in  a  weather  central. 

It  was  further  concluded  that  the  time  required  for  the  prepa- 
ration of  the  complete  forecast  was  within  the  capabilities 
of  the  average  weather  detachment* 

Because  of  its  assumptions  of  barotropy,  the  Pjortoft 
method  is  incapable  of  forecasting  the  changes  in  the  in- 
tensities of  circulation;  hence  it  will  not  show  the  rapid 
and  spectacular  deepening  of  lows  aloft,  which  are  ascribed 
to  the  baroclinic  state  of  the  atmosphere.   This  paper  is  an 
attempt  to  develop  statistical  methods  which  will  indicate 
the  magnitude  of  baroclinic  changes  in  the  upper  air  and, 
together  with  the  Pjortoft  procedures,  provide  an  improved 


objective  forecast  routine  for  the  meteorologist  in  the 
field. 


CHAPTER  II 
THEORETICAL  DEVELOPMENT 
The  vorticity  equation  Is  expressed  in  isobarie 
coordinates  as  follows: 

where  <y  is  the  vertical  component  of  the  relative  vorticity, 
f  is  the  ooriolis  parameter,  V  is  the  horizontal  velocity, 
p  is  the  pressure,  and  CO  -  Jgj*  .  Vectors  are  indicated  by 
overlining,  and  the  vector  operator  V  =  j*  *     *  5i  1   ~*  ^  ^ 
has  its  usual  vector  meaning*  Expanding  the  left  hand  side 
of  this  equation, 

f 
and  combining  with  the  right  hand  side  we  have: 

or  the  local  rate  of  change  of  relative  vorticity  is  equal 
to  the  sum  of  terms  representing  horizontal  adveotion  of 
relative  vorticity,  horizontal  advection  of  the  change  of 
ooriolis  parameter  with  latitude,  vertical  advection  of 
absolute  vorticity,  generation  of  absolute  vorticity  due  to 
divergence  in  the  velocity  field,  and  generation  of  absolute 
vorticity  by  the  wtiltingM  term*   If  we  let  the  sum  of  these 


last  four  terms  be  equal  to  R,  this  equation  may  be  ex- 
pressed as 


J£---V-t/**  '  (4) 


FJortoft,  In  his  simplest  model,  considered  R  negligible 
and  forecast  the  local  rate  of  change  of  relative  vorticity 
by  advecting  the  relative  vorticity  at  the  mean  geos trophic 
wind  velocity.   If  statistical  methods  were  utilized  to 
determine  the  quantity  R  and  this  value  added  to  the  adveoted 
field,  the  forecast  would  be  improved. 

The  purpose  then  is  to  develop  an  objective  prognostic 
routine  consisting  of  the  following  steps:  1)  determine  the 
relative  vorticity  field;  2)  advect  this  field  with  the  mean 
geostrophic  wind;  3)  graphically  calculate  the  change  in 
relative  vorticity  due  to  adveotion;  4)  add  to  this  change 
field  an  empirically  determined  correction  to  compensate  for 
'  the  neglected  terms  in  the  vorticity  equation;  and  5)  recover 
the  height-change  field  and  the  contour  field  from  the 
corrected  vorticity -change  field* 

Well  known  methods  j_5  J  are  available  for  graphically 
computing  the  field  of  relative  vorticity.  Making  use  of 
the  geostrophic  approximation,  the  relative  vorticity  can 
be  expressed  as 

where  Z  is  the  contour  height  and  *\J     is  the  Laplaolan 
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operator  in  the  plane.   Written  in  finite  difference  form 

we  have  a  means  for  determining  the  field  graphically,  where 
the  subscripts  denote  values  of  Z  at  the  four  sides  and  at 
the  center  of  a  square  grid,  and  d  is  the  grid  distance.   If 
the  small  variation  of  f  with  latitude  is  ignored,  the  rela- 
tive vorticity  at  any  point  may  be  represented  by 


where 


J  ■-  *?  - 


*  S   ljr[*.     ***     **J     <*?y] 


(7) 


is  the  space -me an  contour  height. 

Now  combining  equations  (4)  and  (7),  and  replacing  V 
by  the  more  conservative  velocity  field  V*     ,   we  have 

Ar(2-z)    =-$$  .  V„(5-e)  +K  (8) 

where  V&     is  the  geostrophic  mean  wind. 

If  the  equation  for  relative  vorticity  (5)  be  differen- 
tiated with  respect  to  time,  we  have 

y.-J>7^.  (9) 

The  solution  to  this  partial  differential  equation  offers 
the  means  of  recovering  the  height-change  field  from  the 
forecast  vorticity-change  field.  PJortoft  TlJ  has  developed 
a  graphical  solution  to  this  equation  by  expanding  in  a  con- 
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verging  infinite  series  and  then  neglecting  the  smaller 
terms.   If  all  but  the  first  two  terms  are  neglected,  the 
height  change  is  approximated  by 

where  ^   denotes  the  24-hour  change  and  the  bar  denotes  the 
space-averaged  value  of  the  quantity.  For  a  less  accurate 
approximation,  all  but  the  first  term  may  be  ignored  giving 

4  2  --  -A/  (Jj) 

or  the  24-hour  height  change  is  equal  to  the  negative  of  the 
24 -hour  vorticity  change. 


CHAPTER  III 
PROCEDURE 
Used  in  this  research  was  a  series  of  500-mb  charts  of 
the  United  States  and  the  Pacific  Ocean  prepared  by  the  staff 
analysts  at  the  U.  S.  Naval  Postgraduate  School,  Monterey, 
California,  during  the  months  of  February  and  March  1957. 

The  field  of  24-hour  vorticity  change  due  to  advection 
was  determined  by  the  Pj or toft  method.   This  field  was  then 

s 

subtracted  from  the  actual  vorticity-change  field  as  calcu- 
lated from  the  next  day*s  actual  vorticity.   This  field,  R 
in  equation  (8),  was  considered  to  be  the  vorticity  change 
due  to  baroclinic  influences  in  the  atmosphere.   The  fore- 
casting of  this  quantity  will  be  considered  in  the  next 
chapter. 

In  constructing  vorticity  charts  by  graphical  means 
a  grid  distance  of  6  degrees  of  latitude  at  60  degrees  North 
was  considered  to  give  best  results.   Isopleths  of  relative 
vorticity  were  advected  at  100  percent  of  the  space-mean 
500-mb  wind.  Most  authors  (~4,  6l  suggest  80  percent  as 
being  best  in  order  to  compensate  for  the  use  of  the  500-mb 
level  as  the  level  of  non- divergence.   However,  as  the  objective 
of  this  paper  is  to  forecast  the  effects  of  divergence  on  the 
vorticity  field,  100  percent  of  the  mean  wind  was  used. 

Some  difficulty  was  encountered  in  the  graphical  tech- 
niques of  adding  and  subtracting  due  to  fields  consisting  of 
many  small,  closed  centers.   This  difficulty  was  usually 
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overcome  by  drawing  intermediate  iaopleths.   However,  on  a 
few  occasions  the  field  was  still  not  adequately  determined; 
but  it  was  evident  that  this  occurred  only  when  the  vorticity 
change  was  small  and  due  entirely  to  advection. 
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CHAPTER  IV 
FORECASTING  VORTICITY  CHANGES 

As  discussed  in  Chapter  II,  two  of  the  sources  of  non- 
adveotive  relative  vorticity  changes  are  divergence  in  the 
velocity  field  and  vertical  motions  in  the  atmosphere.   An 
attempt  will  be  made  to  forecast  these  changes  by  purely 
statistical  means, 

Scherhag  ^7  J  states  that  where  upper  contours  converge, 
an  "entrance  zone",  pressure  rises  will  occur  due  to  trans- 
verse convergence,  and  where  upper  contours  diverge,  a  "delta 
zone",  pressure  falls  will  occur  due  to  transverse  divergence. 
O'Connor  [_8J  states  that  convergence  and  upper  height  rises 
are  associated  with  low  speed  winds  approaching  straight  or 
cyclonic ally  curved  strong  contour  gradients  or  with  high 
speed  winds  approaching  anticyc Ionic ally  curved  weak  contour 
gradients,  and  that  divergence  and  upper  height  falls  are 
associated  with  high  speed  winds  approaching  oyclonically 
curved  weak  contour  gradients.   These  two  statements  provide 
the  parameters  that  will  be  correlated  with  non-adveotive 
vorticity  changes  in  order  to  forecast  the  magnitude  of  the 
changes  due  to  divergence  in  the  velocity  field. 

O'Neill  [_9  J  conducted  an  investigation  as  to  the  validity 
of  the  above  statements  for  forecasting  height-change  centers 
and  concluded  that,  although  qualitatively  correct,  the  rules 
devised  did  not  show  sufficient  correlation  between  observed 
and  forecast  height  change  to  provide  a  quantitative  fore- 
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cast.   However,  It  is  felt  that  even  though  no  high  corre- 
lation existed  between  total  height-change  and  selected 
parameters,  such  a  correlation  may  exist  when  non-advectlve 
vortlclty  change  Is  the  dependent  variable. 

An  attempt  was  made  to  measure  the  radius  of  curvature 
as  one  of  the  parameters  to  be  used.   Because  of  difficulties 
in  determining  where  to  make  the  measurement  this  was  abandon- 
ed in  favor  of  using  five  subjective  classifications:  marked 
cyclonic  curvature,  moderate  cyclonic  curvature,  negligible 
curvature  or  straight  contours,  moderate  anticyclonic  curva- 
ture, and  marked  anticyclonic  curvature.   These  five  cases 
will  be  examined  in  detail. 

Case  I  -  Marked  cyclonic  curvature.   There  were  a  total 
of  21  configurations  that  were  considered  to  fall  within 
this  classification  in  the  series  of  maps  examined.   All 
were  associated  with  major  troughs,  one  located  over  the 
west  coast  of  the  United  States  moving  slowly  eastward  and 
the  other  over  the  Western  Atlantic  which  moved  off  the  chart. 

The  parameter  used  to  measure  the  divergence  of  contours 
is  defined  as  x±     -  d^  -  dj_  where  d^  is  the  distance  in 
degrees  latitude  between  400  feet  of  contour  at  the  trough 
line,  and  d^  is  the  distance  between  400  feet  of  contour  at 
the  upstream  inflection  point  of  the  same  contours.   The 
contours  selected  are  those  in  the  geometric  center  of  the 
flow  around  the  trough.  Figure  1  shows  the  configurations 
involved  and  the  method  of  measuring  x^. 
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Figure   1 
Contour   Pattern    for  Case  I     and 
Method    of  Measuring    -y, 
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Multiple  correlation  techniques  of  R  with  two  independent 
variables,  magnitude  of  the  vorticity  change  and  location 
relative  to  the  trough  line  of  the  vorticity -change  center, 
were  attempted,  but  no  significant  correlation  was  found  to 
exist^.   In  19  cases  there  were  closed  centers  of  non-advectlye 
vorticity  change,  the  mean  position  of  which  was  located 
approximately  on  the  middle  one  of  the  three  involved  contours 
at  a  distance  downstream  from  the  trough  line  of  five  degrees 
of  latitude.   Therefore,  this  was  taken  as  the  position  at 
which  R  was  measured. 

Figure  2  is  the  scatter  diagram  of  the  values  obtained 
in  the  investigation  of  Case  I  examples  with  the  regression 
line  of  R  on  x^,  R  a  116  -f  46  x^.  The  correlation  coefficient 
between  the  two  variables  was  calculated  to  be  0.60. 

It  is  to  be  noted  that  there  were  only  4  cases  in  which 
the  value  of  x^  was  negative,  that  is,  the  geostrophio  wind 
at  the  trough  was  greater  than  at  the  inflection  point.   If 
the  four  cases  were  eliminated,  the  correlation  coefficient 
would  be  improved. 

Case  II  -  Moderate  cyclonic  curvature.  Fifteen  examples 
were  observed  in  this  classification.   These  were  areas  of 
large  longitudinal  extent  with  generally  smooth  flow  and 
little  change  in  contour  spacing. 

These  techniques  were  applied  in  all  cases  with  similar 
results. 


is 
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The  values  of  the  parameter  selected  for  the  measurement 
of  divergence,  x2,  were  obtained  as  in  Case  I.   There  were 
only  four  closed  centers  of  non-advective  vorticity  change 
associated  with  this  configuration  of  contours,  the  mean 
position  of  which  was  different  from  that  in  Case  I;  however, 
in  the  interest  of  consistency,  the  same  location  for  measure- 
ing  R  was  used  in  both  cyclonic  cases. 

Figure  3  is  the  scatter  diagram  for  the  15  Case  II  ex- 
amples,  A  regression  line  of  R  on  x^f   R  a  24  -f-  78  x2#  was 
obtained.   These  data  gave  a  correlation  coefficient  of  0.71. 

Case  III  -  Negligible  curvature.   There  were  18  examples 
observed  under  this  classification.   They  were  looated  in 
broad  troughs  or  ridges  and  in  areas  of  extensive  zonal  flow. 

The  value  of  the  parameter  X3  is  defined  as  x3  s  d2  -  d-^ 
where  d-^  is  the  distance  in  degrees  latitude  between  400  feet 
of  contour  at  the  upstream  end  of  the  straight  contours,  and 
d2  is  the  distance  in  degrees  latitude  between  the  same  con- 
tours at  the  downstream  end.   The  contours  used  were  those  in 
the  center  of  the  flow.   The  value  of  R  was  measured  at  the, 
downstream  end  of  the  center  contour.  Figure  4  indicates  the 
parameters  used. 

The  scatter  diagram,  figure  5,  gives  the  results  of  the 
Case  III  examples  with  the  regression  line  of  R  on  x$, 
R  s   19  +  45  X3,  with  a  correlation  coefficient  of  0.79. 

Case  IV  -  Moderate  anticyclonic  curvature.   There,  were 

15  examples  in  this  category,  looated  in  areas  of  broad  zonal 
flow.   There  was  little  change  in  contour  spacing  involved 
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In  these  areas* 

The  value  of  the  divergence  parameter,  X4,  is  defined 
as  X4  ■  dp  -  dj_,  where  dp  is  the  distance  in  degrees  latitude 
between  400  feet  of  contour  at  the  ridge  line,  and  d^  is  the 
distance  in  degrees  latitude  between  the  same  400  feet  of 
contour  at  the  inflection  point  upstream.   The  value  of  R 
was  taken  five  degrees  of  latitude  downstream  from  the  ridge 
line  on  the  contour  of  least  height  value  of  the  three  used. 
Figure  6  indicates  the  method  of  measuring  the  parameters. 

The  scatter  diagram,  figure  7,  indicates  the  results  of 
the  Case  IV  investigation  with  the  regression  line' of 
R  on  X4,  R  ■  -40  +  53x4.   A  correlation  coefficient  of  0.56 
was  obtained. 

Case  V  -  Marked  anticyclonic  curvature.   There  were  24 
configurations  that  were  classified  in  this  category,  all 
directly  upstream  or  downstream  from  the  major  troughs  of 
Case  I. 

The  value  of  the  parameter  X5  was  obtained  in  the  same 
manner  as  X4  in  Case  IV.   There  were  17  closed  centers  of 
non-advective  vorticity  change  involved  with  a  mean  position 
5  degrees  of  latitude  downstream  from  the  ridge  line  near  the 
contour  of  least  height  value.  R  was  measured  at  this  point. 

Figure  8  is  the  scatter  diagram  for  Case  V  with  the  re- 
gression line  of  R  on  X5,  R  •  -155  •+-  6x5.   A  correlation 
coefficient  of  0.10  was  obtained. 

The  results  of  the  above  part  of  the  Investigation  are 
summarized  in  Table  I. 
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Case  I 

Number  of  Examples  21 

Parameter  (°Lat. )  X^  ■  d^  -  d^ 

Regression  Line  R  =  116 +-46x3^ 

Correlation  Coefficient  .60 

Mean  value  of  R  (feet)  165 

Standard  deviation  of  R  (feet)  153 

Case  II 

Number  of  Examples  15 

Parameter  (°Lat.)  xg  ■  d£  -  d^ 

Regression  Line  R  *  24 + 78x2 

Correlation  Coefficient  .71 

Mean  value  of  R  (feet)  20 

Standard  deviation  of  R  (feet)  95 

Case  III 

Number  of  Examples  18 

Parameter  (°Lat.)  X3  -  dg  -  d^ 

Regression  Line  R  ■  19+45x3 
Correlation  Coefficient  .79 

Mean  value  of  R  (feet)  30 

Standard  deviation  of  R  (feet)  50 

Case  IV 

Number  of  Examples  15 

Parameter  (°Lat.)  X4  •  dp  -  d^ 

Regression  Line  R  *  -40  -+-  53x4 
Correlation  Coefficient  .56 

Mean  value  of  R  (feet)  -30 

Standard  deviation  of  R  (feet)  80 

Case  V 

Number  of  Examples  24 

Parameter  (°Lat. )  X5  ■  dr  -  di 

Regression  Line  R  ■  -155  -f  6x5 
Correlation  Coefficient  .10 

Mean  value  of  R  (feet)  -150 

Standard  deviation  of  R  (feet)  130 


TABLE  I 
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As  indicated  earlier,  in  addition  to  the  divergence  in 
the  velocity  field,  there  is  another  source  of  change  in  the 
local  relative  vorticity  that  would  appear  to  lend  itself  to 
statistical  methods:   that  due  to  vertical  motions  in  the 
atmosphere.   An  attempt  was  made  to  derive  a  correction  for 
the  effects  of  vertical  motions. 

Case  I  samples  were  divided  into  two  categories:  Those 
in  which  the  major  trough  was  on  the  lee  side  of  the  conti- 
nental divide  and  those  in  which  the  major  trough  was  on  the 
windward  side.  Lines  of  regression  and  correlation  coefficients 
of  R  were  computed  for  each  of  these  two  categories.   There 
were  no  significant  differences.  However,  the  number  in  each 
case  was  very  small.  Results  would  be  more  conclusive  if 
examples  were  taken  of  the  two  categories  in  Case  III,  the 
straight  flow  case,  due  to  its  higher  correlation  factor, 
but  the  series  of  maps  used  in  the  research  was  not  adaptable 
to  such  analysis.   Time  considerations  prohibited  further 
testing. 

Conclusions.   The  correlation  coefficients  in  Cases  I, 
II,  and  III  are  significant  at  the  95%   level.   It  would  be 
indicated  that  in  the  regions  of  cyclonic ally  curved  or 
straight  contours  a  corrective  field  can  be  added  to  the 
advected  field  of  relative  vorticity  that  will  improve  the 
forecast  for  a  24-hour  period.  The  correlation  coefficient 
in  Case  IV  is  significant  at  the  90$  level,  and  in  Case  V 
there  was  no  statistical  significance.  In  the  cases  of  anti- 
cyolonically  curved  contours,  the  radius  of  curvature  of  the 
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contour  appears  to  be  much  more  critical,  but  with  moderate 
curvature  a  correction  can  be  applied  that  will  improve  the 
forecast.  With  marked  curvature  a  parameter  that  more  acou« 
rately  determines  the  deviation  from  the  gradient  wind  will 
have  to  be  used. 

It  should  be  noted  that  the  local  rate  of  change  of 
relative  vorticity  in  the  absence  of  marked  curvature  of 
contours  is  due  primarily  to  advection.   Table  I  shows  the 
small  values  of  the  means  and  standard  deviations  of  non- 
advective  changes  in  these  cases* 
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CHAPTER  V 
INDEPENDENT  TEST  RESULTS 

The  500- mb  prognostic  techniques  developed  in  Chapter  IV 
were  tested  with  a  series  of  five  independent  forecasts  pre- 
pared during  the  month  of  April  1957. 

The  24-hour  relative  vorticity  change  due  to  advection 
was  determined  by  the  Fjortoft  method  using  100$  of  the  mean 
500-mb  wind  as  the  advective  field.   The  24-hour  relative 
vorticity  change  due  to  other  sources  was  prepared  by  plotting 
the  values  of  R  determined  from  the  regression  formulas,  and 
then  subjectively  drawing  the  R  field  over  the  entire  map. 
The  sum  of  these  two  fields  was  considered  the  total  24-hour 
relative  vorticity-change  field,  from  which,  the  24-hour 
height-change  field  was  recovered  by  the  application  of  equa- 
tion (10). 

The  prognostic  charts  prepared  were  verified  on  the  basis 
of  forecast  and  observed  24-hour  height  changes.  Thirty-six 
grid  points  on  each  map,  for  a  total  of  180  oases,  were  used 
in  this  phase.  As  a  basis  of  comparison  of  the  value  of  the 
method  of  prognosis,  forecasts  were  also  prepared  by  the  un- 
modified Pjortoft  method. 

The  root-mean-square  errors  for  the  two  methods  showed 
no  significant  difference  when  all  the  180  points  were  con- 
sidered, being  146  feet  for  the  Pjortoft  method  and  142  feet 
for  the  other.  However,  if  only  those  points  were  considered 
for  which  a  value  of  R  had  been  actually  calculated,  a  signifi- 
cant improvement  was  made.   There  were  a  total  of  37  such 
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points  for  the  five  maps.   The  root-mean-square  error  in 
forecast  height-change  for  these  points  was  135  feet  as  com- 
pared to  an  error  of  169  feet  for  the  forecast  prepared  by 
the  unmodified  Fjortoft  method. 

Conclusions.   The  results  of  the  test  forecasts  indicate 
that  the  objective  prognosis  of  the  500-mb  constant  pressure 
chart  cannot  be  accomplished  by  the  methods  developed  herein 
for  the  following  reasons:  1)  subjectivity  enters  into  the 
classification  of  the  contour  pattern;  2)  the  rules  are  appli« 
cable  at  only  a  few  points;  3)  it  is  impossible  to  draw  an 
extended  area  chart  of  the  values  of  R  given  the  values  at 
only  these  few  points;  and  4)  the  method  fails  in  areas  of 
marked  anticyclonic  flow.  However,  for  limited  application, 
it  is  believed  that  the  values  determined  from  the  configu- 
ration of  the  contour  field  can  be  subjectively  employed  to 
improve  the  Fjortoft  prognostic  chart. 
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CHAPTER  VI 
FURTHER  STUDIES 

A  further  study  was  conducted  in  order  to  attempt  to 
discover  continuously  variable  parameters  for  forecasting 
the  non-advective  relative  vorticity  change.  Used  in  this 
phase  of  the  research  was  a  series  of  maps  for  November  1955. 

There  was  no  indication  of  any  continuity  in  the  paths 
of  the  centers  of  values  of  R.  Hence  extrapolation  was  of 
no  prognostic  value. 

A  measure  of  the  baroclinity  of  the  atmosphere  is  the 
number  of  intersections  of  isotherms  and  contours  on  a  con- 
stant  pressure  surface.  An  attempt  was  made  to  use  this 
principle  for  forecasting  in  two  ways.  First,  the  actual 
density  of  the  solenoids  was  plotted  against  the  values  of  R. 
This  was  done  in  several  different  ways,  for  example,  by 
counting  the  number  of  intersections  of  isotherms  and  contours 
within  a  radius  of  three  degrees  of  latitude  of  the  point  at 
which  R  was  measured.   Second,  multiple  correlation  was 
attempted  between  the  gradients  of  temperature  and  height 
value  and  R.   Neither  of  these  attempts  met  with  success. 

Persistence  of  the  values  of  R  was  determined  to  be  the 
most  successful  prognostic  tool  in  the  forecasting  of  R. 
That  is  to  say,  for  forecasting  the  value  of  R  at  any  point 
on  the  chart,  best  results  were  obtained  by  making  the 
assumption  that    ^  r  Qt 
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For  the  seven  maps  considered,  using  100$  of  the  mean 
flow  for  the  advection  of  relative  vorticity,  the  Pjortoft 
method  gave  a  root-mean-square  error  of  211  feet  in  forecast 
height  change  as  approximated  by  equation  (11).  Using  per- 
sistence as  the  forecast  of  R  this  error  was  reduced  to  179 
feet. 

As  mentioned  in  Chapter  I,  Pjortoft  [l]  indicated  a  method 
of  compensating  for  the  change  of  coriolis  force  and  map  scale 
with  latitude.  In  this  improved  model,  equation  (8)  was  re- 
placed by 

where  tf    u  .%.    .      ^  / 

and  sz    is  the  earth1 s  rotational  velocity,  d  is  the  grid 
distance,  m  is  the  map  scale  constant,  (f    is  the  latitude, 
g  is  gravity,  and 

is  a  steady  easterly  wind.  Also  mentioned  earlier  was  the 
evidence  that  80$  of  the  advective  field  gives  better  results. 

The  same  forecasts  were  prepared  using  this  more  refined 
Pjortoft  method  and  the  recommended  80$  wind.   The  root-mean- 
square  error  in  forecast  height  change  was  203  feet.  Using 
persistence  as  the  forecast  for  R,  this  was  reduced  to  177 
feet. 

It  appeared  that  persistence  of  R  relative  to  the  contour 
field,  rather  than  to  the  actual  geographic  position,  might 
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give  superior  results,  but  evidence  was  inconclusive. 

Conclusions,  Any  conclusions  reached  concerning  the 
above  results  must  be  qualified  ones  because  of  the  meager- 
ness  of  the  tests  on  which  they  are  based.  Although  the 
improvement  in  root-mean-square  error  in  the  forecast  height 
change  is  statistically  significant,  only  one  series  of  maps 
was  studied.  A  different  synoptic  situation  may  conceivably 
produce  entirely  different  results.  Further  testing  will  have 
to  be  conducted  before  objective  application  of  the  rules 
would  be  warranted* 

The  primary  effort  in  this  research  was  directed  at 
developing  forecast  rules  based  on  the  configuration  of  the 
contour  field.  The  methods  developed  failed  in  objective 
application  for  the  reasons  given  in  Chapter  V.  As  a  result 
of  the  further  study,  these  objections  have  been  overcome  by 
assuming  the  persistency  of  R.  Preliminary  results  indicate 
that  this  assumption  is  a  valid  one  and  can  be  objectively 
applied  to  prognosis  of  the  500-mb  contours.  Pinal  con- 
clusions can  be  reached  only  after  further  testing. 
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